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SUPERCOOLED  AND  SUPERHEATED  WATER 


C.  A.  Angell 
Department  of  Chemistry 
Purdue  University 
West  Lafayette,  Indiana  47907 


ABSTRACT 

Various  thermodynamic,  transport  and  structural  features  of  vater  in  its 
metastable  states,  both  supercooled  and  superheated,  have  been  determined  in 
recent  years,  and  highly  anomalous  behavior  has  been  found  near  the  limiting 
temperatures  in  each  case.  For  example,  the  heat  capacities  increase  ex¬ 
ponentially  at  each  extreme.  In  the  high  temperature  extreme  this  behavior 
is  predictable  from  van  der  Waals  theory,  but  the  low  temperature  anomaly 
la  not.  The  measurements  will  be  reviewed,  with  emphasis  on  the  more 
extensively  studied  supercooled  states  at  normal  and  high  pressure,  and 
evidence  will  be  given  that  the  low  temperature  anomalies  can  be  described 
by  expressions  containing  a  singular  temperature  closely  associated  with  the 
temperature  of  homogeneous  nucleatlon  of  ice  I. 

INTRODUCTION 

In  comparison  with  the  long-standing  and  intensive  investigations  of 

the  properties  of  water  in  its  thermodynamically  stable  states,  the  interesting 

extensions  above  the  normal  boiling  point  and  below  the  normal  freezing  point 

have  been  more  or  less  neglected.  Although  Kaemtz^1  *  measured  the  vapor 

(2) 

pressure  of  water  at  -19° C  in  1820  and  Regnault  extended  the  measurements  to 
-32. 8° C  in  1847,  there  has  been  a  surprisingly  low  level  of  interest  in 
extending  and  improving  the  extent  of  knowledge  of  the  physical  properties 


of  water  in  this  region.  Even  lass  effort  has  been  devoted  to  characterising 
water  under  superheated  conditions!  although  the  general  pattern  of  behavior 
can  be  determined  from  measurements  made  at  pressures  above  the  saturated  vapor 
pressure  using  short  extrapolations  to  lower  pressures. 

Much  more  effort  has  been  devoted  to  determining  the  limits  of  supercooling 
and  superheating,  particularly  the  former  which  is  of  special  Interest  in 
relation  to  the  Initiation  of  raindrops  by  destabilization  of  clouds  of  super¬ 
cooled  water  droplets.  Thus  It  has  been  known  since  the  time  of  Dufour,  1861, 
that  small  drops  of  clean  water  Isolated  from  solid  surfaces  by  suspension 

between  layers  of  Immiscible  liquids,  can  be  supercooled  consistently  to  at 
(3) 

least  -34°C.  More  recently  meteorologists,  producing  myriads  of  small,  ( ^  1  imi) 
absolutely  clean,  water  droplets  in  cloud  chambers  have  determined  that  there 
Is  a  rather  sharply  defined  limiting  range  of  -41  to  -42°C  in  which  VJOZ  of 
droplets  in  a  given  sample  will  freeze  within  a  period  of  seconds* 

It  Is  assumed  that  this  limit  is  determined  by  a  process  of  spontaneous  order 
fluctuations  intrinsic  to  the  pure  liquid  phase,  as  distinct  from  the  hetero¬ 
geneous  process,  initiated  on  contaminant  solid  surfaces,  by  which  liquids 
normally  crystallize.  It  therefore  sets  the  limit  on  the  temperature  range 
over  which  we  can  hope  to  determine  the  physical  properties  of  liquid  water, 
at  least  by  measurements  conducted  on  normal  time  scales. 

In  the  superheated  region,  in  which  there  is  considerable  technological 
Interest  because  of  explosion  hazards  in  power  plant  heat  exchangers,  determina¬ 
tions  of  limiting  temperatures  have  been  made  by  a  number  of  workers .  ^  These 

have  usually  been  based  on  observations  of  spontaneous  rupture  of  small  water 
droplets  rising  in  liquid  columns  through  a  temperature  gradient.  In  the 
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■ost  successful  studies,  vster  droplets  of  'V'IOvob  diam.  have  been  observed  to 

reach  a  temperature  of  279. 5°C  at  1  atm  pressure  before  explosive  vaporization 
(9,10) 

occurs . 

In  contrast  with  the  case  of  supercooling  the  limit  on  superheating  can 
be  approximately  predicted  from  equations  of  state,  which  imply  the  existence  of 
regions  of  mechanical  instability  in  which  neither  superheated  liquid  nor 
supercooled  vapor  can  exist.  For  instance,  a  van  der  Waals  equation  with 
parameters  based  on  critical  point  data^11^  predicts  the  behavior  shown  by 
the  set  of  isobaric  T-V  plots  in  Fig.  1.  The  maxima  correspond  to  minima 
in  the  corresponding  P-V  isotherms,  hence  to  points  at  which  the  liquid 
becomes  mechanically  unstable  and  must  spontaneously  rupture.  The  maximum 


V  — * 


ig.  1.  V,  T,  isobars  for  water  in  the  van  der  Waals  approximation  showing  vari¬ 
ation  with  pressure  of  spinodal  points  at  limit  of  superheating  of  liquid 
and  supercooling  of  vapor. 


predicted  for  the  1  eta  pressure  Isobar  is  277*C  very  close  to  the  observed 
superheating  limit.  Of  course  the  van  der  Waals  equation  Is  Inadequate  to 
describe  water  and  the  observed  limit  must  correspond  to  the  klnetlcally 
controlled  reflection  of  the  true  stability  limit  which  would  be  located  a 
little  higher  in  temperature.  In  this  region  the  Instability  must  be  antici¬ 
pated  by  large  fluctuations  in  density  and  entropy  and  these  will  be  reflected 
in  the  magnitudes  of  the  related  thermodynamic  properties,  compressibility  and 

heat  capacity,  which  we  will  discuss  below. 

There  is  no  corresponding  accepted  body  of  theory  which  predicts  a 
mechanical  stability  limit  for  the  supercooling  temperature  range  for  liquids, 
but  the  correspondences  we  will  describe  between  liquid  property  anomalies  in 
the  two  extreme  ranges  will  suggest  that  the  development  of  such  a  theory 
may  be  needed,  at  least  for  the  case  of  water. 

MEASUREMENT  OF  PHYSICAL  PROPERTIES  IN  METASTABLE  REGIONS 

One  of  the  reasons  for  the  relative  lack  of  information  on  metastable 
region  properties  is  clearly  the  difficulty  of  maintaining  the  sample  under 
study  in  the  (thermodynamically  unstable)  liquid  state  during  any  measurement. 
Obtaining  data  under  these  conditions  requires  the  development  of  special 
techniques.  As  implied  in  the  previous  section,  the  transition  to  the  stable 
state  can  be  impeded  by  removing  all  sources  of  external  surfaces  which 
catalyze  the  transition  to  the  stable  state.  For  macroscopic  samples  this 
requires  exhaustive  cleaning  procedures.  For  purely  statistical  reasons, 
however,  much  less  demanding  sample  preparation  procedures  may  suffice  if 
very  small  samples  are  used.  Use  of  small  samples  has  been  the  principle  on 
which  most  of  the  recent  measurements  of  water  properties  under  metastable 
conditions  have  been  carried  out. 
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There  are  three  approaches  to  small  sample  preparation.  One  may  use 

(a)  very  thin  films  (between  plastic  or  glass  sheets),  for  e.g.,  spectroscopic 

studies  in  strongly  absorbing  regions; (b)  small  diameter  columns  (in 

glass  or  plastic  capillaries)  for  volumetric,  diffusion,  or  ESR  studies; ^  ^ 

(c)  very  short  small  diameter  columns  or  droplets,  e.g.  for  NMR,  viscosity, 

(19-23) 

heat  capacity  and  some  spectroscopic  studies.  For  each  dimension 

made  microscopic  in  the  above  series,  a  gain  in  limiting  temperature  is  made 
for  a  sample  of  given  intrinsic  cleanliness.  With  droplets,  short  term  measure¬ 
ments  can  be  made  down  to  -38°C  at  the  low  temperature  and  up  to  ^  270°C  in 

the  high  temperature  region.  Many  important  results  have  recently  been  obtained 

(22) 

using  the  Rasmussen  MacKenzie  emulsification  technique  in  which  water  is 
dispersed  in  'V  3  ym  droplets  in  heptane  or  (similar  hydrocarbon)  supersaturated 
with  SPAN  65  (sorbitar.  tristearate). 

The  actual  methods  in  which  the  small  samples  are  utilized  in  each 
particular  measurement  will  not  be  reviewed  here  for  lack  of  space,  and  the 
reader  is  referred  to  the  original  literature  cited  in  the  following  sections. 

PHYSICAL  PROPERTIES  OF  METASTABLE  WATER 
Superheated  Region 

DIRECT  MEASUREMENT:  SOUND  VELOCITY  AND  ADIABATIC  COMPRESSIBILITY 

It  appears  that  despite  the  number  of  attempts  to  determine  the  super¬ 
heating  limit,  the  only  physical  property  directly  measured  under  superheated 

(24) 

conditions  is  the  sound  velocity.  This  has  been  studied  by  Trinh  and  Apfel 

to  'v  180°C,  which  is,  however,  far  short  of  the  limit  at  'u  280°C.  The  results 

at  the  highest  temperature  differ  slightly  from  those  of  McDade  et  al  obtained 

(25) 

under  saturation  conditions,  but  concur  in  the  continued  curvilinear  decrease 
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These  are  compared  with  the  isothermal  compressibilities  extrapolated 
by  Kell  and  Whalley  in  Table  1  which  is  taken  from  Ref.  24. 


TABLE  1  Thermodynamic  parameters  of  superheated  water  under 
atmospheric  pressure  (after  Trinh  and  Apfel  (1978)). 


T 

(°C) 

P 

(g/cm3) 

c 

(m/sec) 

e  x  io11  eT  x  io11 

8  .  2  * 

(cm  /dyne) 

y 

110 

0.9509 

1,532 

4.45 

5.08 

1.142 

120 

0.9428 

1,518 

4.60 

5.31 

1.154 

130 

0.9344 

1,501 

4.75 

5.57 

1.172 

140 

0.9256 

1,482 

4.92 

4.87 

1.193 

150 

0.9165 

1,460 

5.11 

6.21 

1.215 

160 

0.9070 

1,437 

5.34 

170 

0.8968 

1,412 

5.59 

EXTRAPOLATED  PROPERTIES 

Since  accurate  measurements  of  physical  properties  can  be  made  at  the 
saturation  vapor  pressure  and  at  higher  pressures,  and  since  the  saturation 
vapor  pressure  does  not  exceed  even  10  atm  before  180°C,  many  properties 
of  superheated  water  at  1  atm  pressure  can  be  estimated  quite  accurately  by 
short  extrapolations  vs.  pressure.  Such  extrapolations  have  been  made  to 
150#C  by  Kell  and  Whalley  for  the  specific  volume,  using  data  of  their  own  and  of 
other  authors  as  described  in  their  definitive  paper,  ref.  26.  Data  at  0 
and  10  bar  (0.1  MPa)  applied  pressure  are  reproduced  in  Table  2. 

TABLE  2  Specific  volume  of  superheated  water  to  150*C, 

(Kell  and  Whalley,  1965) 


T  Density 


0.01 

MPa 

0.1 

MPa 

100 

1.043451 

1.042992 

110 

1.051594 

1.051114 

120 

1.060364 

1.059858 

130 

1.069791 

1.069256 

140 

1.079900 

1.079331 

150 

1.090735 

1.090126 

▲  method  for  estimating  the  specific  heet  of  water  for  temperatures 


outside  the  stable  region  has  been  described  by  Llenhard 
the  expression 


(27) 


who  di  _ves 


C  (T  ,p  )  -  C  + 
p  r  r  pref 


0.001644  (1-1. 2p) 


(T  -  T  ) 
v  r  r 

m 


T7T 


(2) 


where  T  ■  T/T  ,  T  the  critical  temperature  (-  647.2  K) 
r  c  c 

P  ■  p/p  »  P  the  critical  pressure  (»  217.7  atm) 
r  c  c 

C  "Is  the  heat  capacity  at  arbitrarily  chosen  reference  pressure, 
pref  e.g.  200  atm 


■  the  maximum  superheating  temperature  at  p^ 


Using  the  ASME  Steam  Table  data  for  C  at  the  reference  pressure  and  taking 

P 

T  at  279. 5°C,  Eq.  (2)  yields  the  behavior  shown  in  Fig.  3.  The  plot  is 

rm 

presumably  in  error  at  the  high  temperature  extreme  because  the  splnodal  at 
1  atm  pressure,  where  C  •►  *  ,  will  exceed  the  measured  superheating  limit 
279. 5°C  by  a  short  but  unknown  interval.  The  full  curves  show  the  predicted 

C  for  the  case  T  -  290  and  299  K. 

P 

In  either  case  the  exponential  Increase  in  ,  as  the  superheating 
limit  of  279. 5°C  is  approached,  is  striking  and  should  be  borne  in  mind 
in  considering  the  observations  on  supercooled  water  in  the  following  section. 

Superheated  water  has  been  observed  recently  at  temperatures  as  low  as 
-18®C  by  Henderson  and  Speedy,  using  a  technique  which  places  small 
samples  of  water  confined  in  capillaries  under  isotropic  tension.  Superheated 
water  at  a  pressure  of  -118  bar  (—1.18  MPa)  has  a  density  maximum  at  6.1°C 
rising  to  6.7#C  at  -156  bar.  At  the  lower  temperatures  of  their  study,  the 
water  sample  was  simultaneously  superheated  and  supercooled. 
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Fig.  3.  Heat  capacity  of  superheated  water  at  one  atm  pressure  from 

Eq.  (2)  with  T  (spinodal)  taken  to  be  either  279.5  (curve  1) , 
rm 

290  (curve  2)  or  300  (curve  3) . 
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Supercooled  Region 

THERMODYNAMIC  PROPERTIES 

The  thermodynamic  property  which  has  been  measured  over  the  widest 
temperature  range  at  normal  pressure,  and  which  now  has  been  measured  with 
perhaps  the  greatest  accuracy,  is  the  constant  pressure  heat  capacity. 

It  la  probably  true,  also  that  the  behavior  of  this  fundamental  property  of 
water  Is  the  most  immediately  surprising  of  the  many  unusual  aspects  of  the 
low  temperature  regime  for  this  liquid,  so  it  is  appropriate  to  present  these 
results  first  in  this  article. 

Though  known  since  1890  to  -5°C  (Martinetti) ,  the  heat  capacity  was  only 

shown  to  be  highly  unusual  in  the  supercooled  state  by  the  emulsion  measurements 
(29) 

of  Rasmussen  et  al  which  extended  to  -38°C.  These  authors  presented  data 
from  both  differential  scanning  calorimetry  and  drift  calorimetry,  but  the 
measurements  were  preliminary  and  not  of  great  accuracy .  There  has  been 
some  effort  recently  in  the  author's  laboratory  to  improve  the  accuracy  and 

precision  of  this  measurement.  Using  the  Perkin  Elmer  DSC-2  differential 
scanning  ca lo rime te r ,  sets  of  data  have  been  obtained  by  different  operators, 

using  different  emulsification,  calibration,  and  water  content-determining 
procedures,  which  are  self-consistent  to  ±2Z.^  ^  These  data  are  presented 

in  Table  3  and  Figure  4.  Some  data  at  higher  pressures  have  been  derived 
by  Kanno  and  Angell^^  from  combination  of  the  Table  3  data,  and  high  pressure 
volumetric  results,  discussed  below,  using  the  thermodynamic  relation, 

V>  •  yo)  -  t  j[  (!p)  <» 

However,  these  are  not  considered  of  sufficient  accuracy  to  present  in  tabular 
form. 
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TABLE  3  Heat  capacity  of  supercooled 


t*C 


water  at  one  atm  pressure 


Specific  Heat  (J  mole  *  deg 


Bulk  Sample 


Emulsion  Sample 
(Average) 


75.7 

76.0 

76.6 

76.2 

77.6 

77.4 

78.9 

81.4 

85.4 

92.7 

96.7 

-20 


T/°C 

20 


40  60  80 


•  •  •  BULK  SAMPLE 
ao°  MCH  CARRIER 

•  •  •  N-OCTANE  CARRIER 

•  •  •  LITERATURE  DATA 


Density  measurements  have  been  performed  for  H^O  and  D^O  using  both 

(32  33) 

capillary  column  length  measurements  (by  Schufle  and  coworkers  ’  and 

(34)  (22) 

Zheleznyi  )  and  emulsion  volume  measurements  (by  Rasmussen  and  Mackenzie  ') 

with  results  in  good  accord.  Specific  volumes  and  derived  expansivities  extend 

down  to  -34°C  at  one  atmosphere  pressure.  These  are  presented  in  Table  4. 


Table  4 

Density,  Molar  Volume  and  Expansivity  of  H^O  and  D^O  at  One  Atm  Pressure 

(after  Zheleznyi,  1969) 


h2° 

d2° 

t*C 

density 
(p/g  cc  1) 

molar  volume 

(V/cc  mol  ) 

expansivity 
10^a/deg  1 

density 

p/g  cc 

molar  volume 
V/cc  mol  * 

expansivity 
10^a/deg  1 

0 

0.9999 

18.018 

1.1030 

18.131 

-  5 

0.9995 

18.025 

-  1.7 

1.1013 

18.160 

-  4.12 

-10 

0.9983 

18.046 

-  3.1 

1.0987 

18.204 

-  5.9 

-15 

0.9964 

18.080 

-  4.7 

1.0950 

18.265 

-  8.5 

-20 

0.9936 

18.130 

-  7.0 

1.0893 

18.360 

-13.0 

-25 

0.9895 

18.206 

-10.8 

1.0809 

18.503 

-30 

0.9829 

18.329 

-16.5 

-34 

0.9751 

18.474 

' 

Using  a  high  pressure  capillary  technique,  these  measurements  have  been  extended 

(31) 

to  190  MPa  by  Kanno  and  Angell,  but  only  over  a  more  limited  temperature 
range  and  only  for  D20.  Their  data  are  shown  in  Table  5  which  includes  data 
from  Kanno  and  Angell  (1979)  for  normal  pressure.  Note  how  the  temperature 
of  the  density  maximum  is  depressed  by  increasing  pressure. 


ISOTHERMAL  COMPRESSIBILITY  K.  /  MPa'1 


Fig.  5.  Isothermal  compressibilities  of  water  at  normal 

and  high  pressures,  from  Kanno  and  Angell  (1979). 


15 


Adiabatic  compraaalbllltlea  may  be  calculated  from  the  above  data  ualng 
tha  relation 


K  C  /C 
T  v  p 


q2VT 

c 

p 


(4) 


or  obtained  directly  from  sound  velocity  measurements.  The  latter  have  been 

obtained  by  several  groups  working  with  different  techniques  in  both  ultrasonic 

and  hypersonic  frequency  ranges  5 36tn?each  frequency  range  there  ia  essential 

agreement  amongst  different  groups  down  to  -15#C,  but  there  is  a  persistent 

and  strange  difference  between  the  ultrasonic  data  on  the  one  hand,  and  the 

(38) 

aero  frequency  (Eq.  4)  and  the  hypersonic  data  on  the  other.  Leyendekkers 
has  reviewed  the  area  and  concludes  that  the  isoentropic  assumption  central  to 
the  data  analysis  is  inadequate  in  the  case  of  ultrasonic  pressure  waves. 

Data  are  presented  graphically  in  Fig.  6.  Below  -15#C  there  are  real  data 
discrepancies  amongst  different  workers  which  remain  to  be  resolved. 


Fig.  6.  Adiabatic  compressibilities  calculated  from  ultrasonic 
sound  velocity  measurements  and  from  Eq.  (4).  At  0°C 
the  latter  value  is  accurately  known,  and  the  large 
uncertainties  indicated  by  the  error  bars  do  not  apply. 
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Vapor  pressures  in  the  supercooled  states  of  water  were  first  measured 

by  Regnault  In  1847  whose  data  extended  to  -32.8*C.  Due  to  Instrumental 

limitations  these  extraordinary  measurements  were  not  of  great  accuracy  and 

(39) 

are  not  quoted  here.  Scheel  and  Heuae  in  1909  obtained  data  to  -15*C 

(40) 

which  have  been  confirmed  recently  by  Bottomley.  A  Van't  Hoff  plot  of 

the  data  between  0  and  -15°C  is  linear,  and  a  vapor  pressure  of  1  mm  Hg  (13.3  Pa) 

la  Indicated  at  -20*C.  Some  data  are  given  in  Table  7. 


Table  7 

Vapor  Pressure  of  Supercooled  Water 
(Data  selection  from  Scheel.  and  Heuse  (1909)) 


t/#C 

p/mm 

-  2.814 

3.724 

-  2.816 

3.726 

-  2.830 

3.726 

-  2.838 

3.724 

-  3.918 

3.432 

-  4,842 

3.209 

-  4.884 

3.202 

-  4.900 

3.193 

-  7.250 

2.675 

-  7.279 

2.669 

-  7.301 

2.665 

-  7.313 

2.661 

-  8.804 

2.373 

-  9.115 

2.315 

-  9.812 

2.191 

-15.176 

1.421 

-15.308 

1.412 
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The  static  dielectric  constant  (relative  permittivity)  of  water  below 

0®C  has  been  measured  by  Hasted  and  Shahidi^  ^  using  an  absolute  method  based 

on  dilute  emulsion  samples  analyzed  using  the  Maxwell-Lewin  formula  for  the 

mean  permittivity  e  of  a  dispersion  of  spheres  of  a  substance  of  dielectric 

m 

constant  in  a  matrix  of  dielectric  constant  .  When  the  volume  fraction 

of  spheres  is  <p>  the  relation  is 


E.  "  3  *  £2 


ei  - 


V2£2  + 


-  e. 


(5) 


Their  results,  shown  in  Fig.  7,  have  been  supported  by  independent  measurements 
(42) 

of  Bodge  and  Angell  based  on  the  Maxwell-Wagner  effect  for  concentrated 
emulsions,  using  an  analysis  which  requires  calibration  with  known  values  at 
some  reference  temperature. 


Fig,  7.  Dielectric  constant  (relative  permittivity)  for 

supercooled  water  (after  Hasted  and  Shahidi,  ref.  41). 
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MASS  TRANSPORT  PROPERTIES 

Tht  Intensive  thermo dyne  ale  properties  discussed  sbove  sre  determined  by 
the  msgnltude  of  fluctustlons  in  extensive  properties  which  sre  chsrscterlstlc  of 
the  liquid  st  equilibrium.  The  trsnsport  properties  likewise  sre  determined  by  the 
growth  and  decay  rates  of  these  same  fluctuations,  hence  are  in  a  sense  also 
equilibrium  properties  of  the  liquid.  However,  they  are  almost  invariably 
determined  by  deliberately  perturbing  the  equilibrium  state  of  the  system 
and  observing  the  rate  of  return  to  equilibrium,  hence  they  are  usually 
considered  separately  as  "non-equilibrium"  properties. 

Fundamentally  related  to  the  process  of  recovery  of  equilibrium  after 
most  types  of  perturbation  is  the  ordinary  (mass)  diffusion  process.  This 
has  been  studied  In  the  supercooled  range  by  ftuppacher^1^  using  the  classical 
radiotracer  dif fusion -out-of-cap illary  technique  and  reaching  a  surprising 
-25*C  and  by  Gillen  et  al^17^  using  the  shorter  time  NMR  spin-echo  technique. 

The  latter  authors  reached  -31°C  with  groups  of  selected  capillary  sample 
containers.  Attempts  to  extend  these  measurements  to  lower  temperatures  using 
pulsed  gradient  method  on  emulsion  samples  have  only  been  partially  successful. 
Data  from  the  preferred  study  of  Gillen  at  al  are  presented  in  Table  8 

after  a  uniform  6.7X  upward  correction  to  normalise  to  the  accepted  diaphragm 

(44) 

method  value  for  bulk  water. 

The  viscosity  of  water  was  first  measured  below  0  C  in  1913  by  White  and 

Twining (45)  whose  data  penetrated  to  -14*C.  Fifty  years  later  Hallett(46)  used 

a  clever  cold  capillary  method  to  extend  the  range  to  -24°C.  Very  recently 

(21) 

Osipov  et  al  have  reached  -35#C  with  a  moving  microscopic  slug  technique 
which  gives  results  in  agreement  to  IX  with  those  of  Hallett  in  their 
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TabU  8 


Self-Diffusion  Coefficients  of  Water 
Into  the  Supercooled  Range 
(After  Gillen,  Douglass  and  Hoch  (1972) 


normalized  to  Mills 

datum  at  25SC) 

t°C 

5  2- 

10  D/cm  sec 

25.1 

2.38 

12.2 

1.69 

2.4 

1.12 

-  9.4 

0.75 

-11.6 

0.675 

-14.4 

0.590 

-17.3 

0.509 

-19.1 

0.467 

-21.3 

0.421 

-23.1 

0.364 

-24.9 

0.343 

-26.8 

0.281 

-28.7 

0.234 

-30.6 

0.200 

eoaon  temperature  range.  With  the  Utter  data, viscosity  becomes  the  moat 
extensively  studied  transport  property  for  water,  though  In  the  case  of  D^O, 
the  deuteron  spin  relaxation  time  data  cover  a  wider  equivalent  range.  The 
viscosity  data  collected  In  Table  9  are  combined  with  the  stable  range  data  in  an 
Arrhenius  plot.  Fig.  8,  to  demonstrate  the  extraordinary  departures  from  Arrhenius 
behavior  exhibited  by  water  In  this  low  temperature  range.  The  apparent  activation 
energy  In  the  loweat  temperature  Interval  has  risen  to  14  kcal/mole  (59  kJ/mole) . 
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Table  9 


Viscosity  of  Water  Below  0*C 


t/'C 

n/cp 

White  &  Twining 

Hallett 

Osipov  and  coworkers 

(191  3) 

(1963) 

(1979) 

0 

1.798 

1.79 

1.80 

-  4.70 

2.121 

-  5 

2.16 

2.14 

-  7.23 

2.341 

-  9.30 

2.549 

-10 

2.66 

2.60 

-15 

3.34 

3.23 

-20 

4.33 

4.36 

-25 

6.45 

-30 

10.2 

-32 

12.7 

t 

-34 

16.4 

-35 

18.7 

The  electrical  conductivity  of  water,  which  proceeds  by  a  proton-hopping 

(Grotthus)  mechanism,  has  not  been  properly  measured  Into  the  supercooled 

region.  The  dielectric  relaxation  time,  which  reflects  the  reorientation 

kinetics  of  the  molecules  and  is  described  extraordinarily  well  by  the  simple 

Debye  relaxation  theory  in  the  stable  range,  has  only  been  measured  to  -8°C, 

where  T  -  2.5  x  10_11  a. ^^According  to  the  Debye  model,  in  which  T 
D  u 

scales  with  n/T  ,  Tp  should  grow  to  1.4  x  10  ^  by  -35 °C. 
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10 


T/«C 


Fig.  8.  Arrhenius  plot  of  viscosity  data  for  normal 
and  supercooled  water. 

Nuclear  spin  relaxation  tinea  can  be  measured  on  very  small  and  quiescent 

pure  water  samples  and,  with  almost  equal  accuracy,  on  capillary  or  emulsion 

droplet  samples,  hence  the  data  for  these  relaxation  processes  are  very 

extensive.  Hindman  and  coworkers have  reported  measurements  down 

to  -31°C  for  170  and  H  relaxation  (where  instrumental  electronic  problems 

2 

were  limiting)  and  down  to  -36°C  for  D  in  D^O.  An  example  of  their  results 
in  which  both  bulk  water  and  emulsion  sample  data  are  included  is  shown  in 
Fig.  9.  For  the  and  most  recently  cases,  our  knowledge  of  nuclear  spin 
relaxation  processes  has  been  extended  most  usefully  to  pressures  of  250  MPa 
(2.5  kbar)  and  temperatures  as  low  as  -86*C  by  Lang  and  Ludemann^20^  using  a 
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1000/ T,  T  IN  DEGREES  K 


Fig.  9.  Arrhenius  plot  of  spin  lattice  relaxation  time  for  ^0  in  normal  and 
supercooled  water  (after  Hindman,  ref.  50). 

high  pressure  strengthened  glass  sample  cell.  Using  e  frequency  of 

a>  ■  6.28  x  10  Hr,  these  workers  were  able  to  observe  a  T,  minimum  (at  -70°C, 

0  1 

200  MPa)  hence  were  able  to  derive,  from  (0  ■  1,  a  reorientation  correlation 

o  0 

-9 

time  T.  of  1.6  x  10  sec,  ^  3  orders  of  magnitude  longer  than  at  25°C.  Water 
under  these  conditions  Is  evidently  a  rather  viscous  liquid,  with  a  viscosity 
of  ^>10  poise  (c.f.  glycerol  under  ambient  conditions,  n  ■  15  poise).  If  the  Debye 
reorientation  time  theory  continued  to  apply  even  approximately.  Reorientation 
times  derived  for  a  range  of  conditions  by  Lang  and  Ludemann  are  shown  In 
Fig.  10.  The  temperature  dependence  of  the  reorientation  time  at  high 
pressure  proves  to  be  distinct  from  that  at  1  atm  pressure  and  Instead  Is 


Reorientation  tines  have  been  Measured  at  one  atmosphere  pressure  for 

dilute  radical  solute  species  in  water  at  temperatures  down  to  -33*C,  by 

ESR  measurements.  Ahn^18^  showed  that  the  reorientation  time  for  ditertiary 

-4 

butyl  nitroxide  at  concentrations  of  10  m  followed  the  viscosity  temperature 
dependence  down  to  -24°C,  and  successfully  predicted  the  then  unknown  viscosity 
of  water  at  -33°C. 

STRUCTURE-RELATED  SPECTROSCOPIC  STUDIES 

Several  spectroscopic  studies  have  been  performed  on  water  in  the  supercooled 

state,  and  others  are  in  progress.  NMR  chemical  shift  measurements  which  average 

overall  configurations,  and  infra-red  measurements,  which  are  "seeing"  a  much 

(23  53) 

shorter  time  scale  structure,  have  been  studied  to  -35°C,  *  but  these  results 

are  not  discussed  here  (see  original  literature  and  review  by  the  author  in 
"Water:  A  Comprehensive  Treatise','  Vol.  7  (Ed.  F.  Franks,  1980)). 

DISCUSSION 

The  results  reviewed  above  offer  a  very  demanding  challenge  to  the 
theoretician . 

The  rapid  variation  in  thermodynamic  properties  at  low  temperatures, 

which  are  so  similar  to  that  of  the  heat  capacity  at  high  temperature  (compare 

Figs.  3  and  4),  suggests  the  possibility  at  low  temperature  of  a  thermodynamic 

singularity  analogous  to  the  instability  (spinodal)  point  at  high  temperatures 

described  in  the  introduction.  This  possibility  is  strengthened  when  the  data 

are  plotted  against  the  familiar  critical  point  temperature  function 

(35) 

e  •  (T  -  T  )/T  where  T  is  the  singular  or  spinodal  point.  This  is 

0  8  8 

shown  for  two  properties,  the  expansivity  a  (a  thermodynamic  property)  and 
the  viscosity  n,  (»  transport  property)  in  Fig.  11.  For  a.  s  smell,  linearly 
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temperature-dependent ,  component  has  first  been  subtracted  in  order  that  the 

component  being  plotted  doea  not  become  zero  at  any  temperature.  This  linear 

component  represents  the  "normal"  (non-crltical  or  non-anomalous)  part  of  the 

thermodynamic  property,  which  is  always  superimposed  on  the  anomalous  part  in 

any  static  property  of  the  system.  Such  a  subtraction  need  not  be  made  in  the 

case  of  the  viscosity  because  the  shear  relaxation  time  la  dominated  by  the 

decay  time  of  the  anomalous  fluctuations.  Both  properties  give  linear  plots 

when  T  is  fixed  at  228  K.  Currently,  attempts  are  being  made  to  derive,  from 
8 

binary  solution  studies  in  which  the  anomalous  behavior  at  low  temperatures  is 
eliminated,  appropriate  background  ("normal")  components  for  the  heat  capacity 
mil  isothermal  compressibility  (Angell  and  Tucker,^  ^  Ogunl  and  Angell^^)  so 
that  the  exponents  describing  the  singular  part  can  be  determined. 

Tt*C 


log  € 

Fig.  11.  Plot  of  anomalous  components  of  selected  supercooled  and 

superheated  water  thermodynamic  properties,  and  of  viscosity 
for  supercooled  water  vs  log  (T/T  -  1)  with  T  “  228  K. 

S  8 
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It  la  to  be  noted  that  228  K  la  only  some  5  K  below  the  lowest  .'emperature 

TH*  to  which  »««11  samples  can  be  supercooled  before  crystallisation  occurs . 

This  suggests  strongly  that  the  singular  point  is  indeed  a  low  teiq>erature 

mechanical  stability  limit  for  the  liquid  phase  analogous  to  that  for  the 

superheated  liquid  discussed  earlier,  although  "internal"  liquid  phase 

interpretations  of  Tfl  are  also  possible. 

It  is  unfortunate  that  no  directly  measured  high  temperature  data 

near  the  high  temperature  are  available  to  include  for  comparison  with 

the  supercooled  water  data.  In  their  absence,  we  use  the  Eq.  (2)  -  based 

data  displayed  in  Fig.  3,  after  subtracting  a  constant  "background"  (normal 

component)  of  75  J  mol  ^  deg  \  to  obtain  the  plot  labelled  C  (superheated) 

P 

in  Fig.  11.  The  data  at  large  £  (where  the  anomalous  component  should  become 
dominant,  and  the  error  in  choice  of  background  be  least  Influential)  are 
linear  in  £  with  a  slope  of  1.5.  This,  however.  Is  larger  than  the  expected 
classical  value  1.0,  and  larger  than  even  the  lattice  gas  value,  1.25,  which 
might  be  expected  to  apply  very  near  the  splnodal  if  measurements  were  possible. 
It  is  probable  that,  close  to  T>t  the  form  of  Iq.  (2)  becomes  Inadequate. 

If  we  take  T^  at  the  high  temperature  limit  to  be  300*C,  then  at  279. 5*C 
where  spontaneous  rupture  is  observed,  the  value  of  £  is  almost  the  same 
as  that  of  supercooled  water  at  -40“C  where  homogeneous  nucleation  of  similar 
size  droplets  is  observed. 

Considering  the  low  temperature  behavior  further,  Kanno  and  Angell^^ 

have  shown  that  the  close  relation  between  TR  and  is  maintained  as  pressure 

la  Increased  to  200  MPa  even  though  T  decreases  by  more  than  50  K 

n 

(to  -92  K  at  200  MPa).  On  the  other  hand,  both  thermodynamic  and  relaxa¬ 
tion  time  data  at  pressure  above  200  MPa  indicate^15^  the  disappearance  of 
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singular  behavior  (Tfl  becomes  too  low)  and  the  return  to  "normal"  low  taapera- 
ture  liquid  behavior  in  which  the  temperature  Tq  appears  as  some  kind  of 
Internal  equilibrium  low  temperature  limit  on  the  liquid;  in  effect,  a 
temperature  at  (or  before)  which  it  would  undergo  a  second  (or  higher)  order 
tranaition  to  an  "ideal  glass"  phase. (56-59)  In  practice  liquids  always  pass, 
for  kinetic  reasons,  through  a  non-equilibrium  transition  (the  1  glass  transition  ) 
at  some  temperature  T  a  little  higher  than  Tq.  (So  far,  not  even  the  glass 
transition  has  been  observed  for  water  because  of  fast  crystallization  even 
In  the  region  above  200  MPa,  though  several  attempts  have  been  made  in  the 
author's  laboratory  to  suppress  the  latter  by  fast  quenching  procedures) . 

The  above  observations  on  the  thermodynamic  characterisation  of  water 
are  summrized  In  Fig.  12  in  the  fora  of  a  qualitative,  but  realiatically 


Fig.  12.  Free  energy  surfaces  for  solid,  liquid  and  gaseous  phases 
of  water  showing  metastable  surfaces  and  their  stability 
edges,  as  described  in  text. 
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scaled,  free  energy  surface  diagram  for  water.  This  shows  the  Intersections 

with  stable  crystal  polymorph  and  gas  free  energy  surfaces  which  give  us  the 

normal  crystallization  and  boiling  phenomena,  but  also 

(a)  the  metastable  extensions  of  the  liquid  and  gas  surfaces  at 

T  <  T  ,  P  <  P  ,  showing  the  surface  edges  corresponding 
c  c 

to  mechanical  stability  limits  on  vapor  supercooling  T  and  liquid 
superheating  T  . ; 

(b)  the  metastable  extensions  of  crystal  and  liquid  surfaces 

showing  a  (necessary  but  unobservable)  superheating  limit  T  „ 

for  the  crystal  due  to  shear  mode  Instabilities  at  an  unknown 

high  temperature,  and  the  putative  mechanical  stability 

limit  on  liquid  supercooling  T  , ; 

s,l 

(c)  the  vanishing,  at  higher  pressures,  of  the  low  T  instability  edge 

as  T  ,  falls  below  T  .  In  the  higher  pressures  the  amorphous 
s,l  o 

phase  surface  is  seen.  Ideally,  as  continuing  to  0  K  with  a 
slope  discontinuity  at  TQ  such  that  below  the  surface  has 
essentially  the  same  temperature  derivative  (i.e.  entropy)  as 
the  stable  crystalline  phase  at  that  temperature. 

Providing  a  theoretical  basis  for  the  impending  low  temperature  catastrophe 
la  a  somewhat  difficult  problem,  particularly  as  there  are  no  other  examples 
of  liquids  exhibiting  this  type  of  behavior  known  at  this  time.  The  facts 


(1)  the  anomalies  are  suppressed  by  constraining  the  system  to  constant 


volume  (Cv  actually  decreases  slightly  with  decreasing  temperature  in  the 

(53) 

temperature  range  0  to  -40*C  ),  and 
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(11)  the  anoaallea  are  banlahad  aven  from  constant  pressure  properties 


at  high  pressures  where  the  water  molecules  sre  forced  Into  closer  packing 
arrangements , 

Imply  that  geometrical  aspects  of  the  packing  of  water  molecules  in  the 
preferred  open  tetrahedral  network  configuration  are  of  central  importance. 

IWo  theoretical  approaches  consistent  with  this  obserwtlon  are  currently 
under  development.  One,  due  to  Stanley  and  Teixeira^*^'^  Is  basically  a 
theory  of  gelation  of  fully  bonded  (hence  open -packed)  water  molecules,  in 
which  a  percolation  threshhold  near  -45 °C  appears  to  arise  as  a  natural 
consequence  of  the  random  formation  of  hydrogen  bonds  of  fixed  energy.  The 
other,  by  Stllllngerv  Is  a  more  general  development  of  the  statistical 
mechanics  of  finite  systems  in  which  the  differential  geometry  of  a  con¬ 
figuration  space  surface  seems  to  lead  naturally  to  a  catastrophe  at  sufficient 
supercooling. 

These  ideas  are  yet  at  an  early  stage  of  their  development  and  the 
appropriate  role  played  in  the  phenomenology  by  cooperative  hydrogen  bonding 
(l.e.  energetic  factors  as  opposed  to  purely  geometrical  factors)  imposing 
a  superstructure  of  consequences  on  a  framework  of  random  hydrogen  bonding 
remains  to  clearly  delineated.  It  is  certain,  however,  that  from  a 
successful  treatment  of  the  supercooled  water  phenomenology  there  will 
emerge  a  much  more  complete  and  sophisticated  picture  of  the  nature  of 
liquid  water  than  Is  currently  available. 
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SUMMARY 


A  review  of  recent  measurements  end  calculations  on  properties  of  water 
In  metastable  superheated  and  supercooled  states  shows  the  existence  of 
Interesting  parallels  In  behavior  which  are  probably  unique  to  pure  water 
at  relatively  low  pressure  200  MPa).  For  this  case  the  singular  behavior 
expected  for  all  liquids  at  large  superheating,  as  they  approach  their 
mechanical  stability  limits ,  is  repeated  In  the  low  temperature  supercooled 
regime.  Available  data  are  consistent  with  a  singular  temperature  of  228  K,  at  1  atm. 
More  accurate  and  more  extensive  data  are  needed  to  elucidate  this  behavior. 

The  singular  behavior  Is  eliminated  by  sufficient  Increases  In  applied 
pressure,  as  the  hydrogen  bond  angle  relationships  essential  to  the  manifesta¬ 
tion  of  the  singularity,  are  suppressed. 
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